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ABSTRACT 

The EFAR project is a study of 736 candidate early-type galaxies in 84 clus- 
ters lying in two regions towards Hercules-Corona Borealis and Perseus-Cetus at 
distances cz pa 6000 — 15000 km/s. In this paper we describe a new method of 
galaxy photometry adopted to derive the photometric parameters of the EFAR 
galaxies. The algorithm fits the circularized surface brightness profiles as the 
sum of two seeing-convolved components, an i? 1 / 4 and an exponential law. This 
approach allows us to fit the large variety of luminosity profiles displayed by the 
EFAR galaxies homogeneously and to derive (for at least a subset of these) bulge 
and disk parameters. Multiple exposures of the same objects are optimally com- 
bined and an optional sky-fitting procedure has been developed to correct for sky 
subtraction errors. Extensive Monte Carlo simulations are analyzed to test the 
performance of the algorithm and estimate the size of random and systematic 
errors. Random errors are small, provided that the global signal-to-noise ratio 
of the fitted profiles is larger than pa 300. Systematic errors can result from 1) 
errors in the sky subtraction, 2) the limited radial extent of the fitted profiles, 
3) the lack of resolution due to seeing convolution and pixel sampling, 4) the 
use of circularized profiles for very flattened objects seen edge-on and 5) a poor 
match of the fitting functions to the object profiles. Large systematic errors are 
generated by the widely used simple R 1 ^ law to fit luminosity profiles when a 
disk component, as small as 20% of the total light, is present. 

The size of the systematic errors cannot be determined from the shape of the 
X 2 function near its minimum because extrapolation is involved. Rather, we must 
estimate them by a set of quality parameters, calibrated against our simulations, 
which take into account the amount of extrapolation involved to derive the total 
magnitudes, the size of the sky correction, the average surface brightness of the 
galaxy relative to the sky, the radial extent of the profile, its signal-to-noise ratio, 
the seeing value and the reduced x 2 °f the fit. We formulate a combined quality 
parameter Q which indicates the expected precision of the fits. Errors in total 
magnitudes Mtot less than 0.05 mag and in half-luminosity radii R e less than 
10% are expected if Q — 1, and less than 0.15 mag and 25% if Q = 2; 89% of 
the EFAR galaxies have fits with Q = 1 or Q = 2. The errors on the combined 
Fundamental Plane quantity FP = \ogR e — 0.3{SB e ), where (SB e ) is the average 
effective surface brightness, are smaller than 0.03 even if Q = 3. Thus systematic 
errors on M TO t and R e only have a marginal effect on the distance estimates 
which involve FP. 

We show that the sequence of R}l n profiles, recently used to fit the luminosity 
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profiles of elliptical galaxies, is equivalent (for n < 8) to a subsample of R 1 ^ and 
exponential profiles, with appropriate scale lengths and disk-to-bulge ratios. This 
suggests that the variety of luminosity profiles shown by early-type galaxies may 
be due to the presence of a disk component. 

Subject headings: galaxies: early type - galaxies: clusters - universe: large scale 
structure - galaxies: peculiar velocities 

1. Introduction 

This is the fourth paper of a series where the results of the EFAR project are presented. 
In Wegner et al. (1996, hereafter Paper I) the galaxy and cluster sample was described, 
together with the related selection functions. Wegner et al. (1997, hereafter Paper II) 
reports the analysis of the spectroscopic data. Saglia et al. (1997, hereafter Paper III) 
derives the photometric parameters of the galaxies. In this paper we describe the fitting 
technique used to derive these last quantities. 

A large number of papers have been dedicated to galaxy photometry. The reader should 
refer to the Third Reference Catalogue of Bright Galaxies (RC3, de Vaucouleurs et al. 1991) 
for a complete review of the subject. By way of introduction we give here only a short sum- 
mary of the methods and tests adopted and performed in the past to derive the photometric 
parameters of galaxies. 

Using photoelectric measurements, photometric parameters have been derived by fitting 
curves of growth. The RC3 values are computed by choosing the optimal curve between a 
set of 15 (for T = — 5 to T = 10, see Buta et al. 1995), one for each type T of galaxies. 
Photoelectric data are practically free from sky subtraction errors (< 0.5%), but can suffer 
from contamination by foreground objects. Typically, 5-10 data points are available per 
galaxy, with apertures which do not exceed 100 arcsec and do not always bracket the half- 
luminosity diameter. Burstein et al. (1987) (who fit the i? 1 / 4 curve of growth to derive 
the photometric parameters of a set of ellipticals) discuss the systematic effects associated 
with these procedures. The total magnitudes Mtot and effective radii R e derived are biased 
depending on the set of data fitted. The errors in both quantities are strongly correlated, 
so that Alogi^ - 0.3A(SB e ) « constant, where (SB e ) = M TOT + 5\ogR e + 2.51og(27r) 
is the average surface brightness inside R e . This constraint (Michard 1979, Kormendy & 
Djorgovski 1989 and references therein) stems from the fact that the product R(I) ' 8 varies 
only by ±5% for all reasonable growth curves (from R l l 4 to exponential laws) in a radius 
range 0.5i? e < R < l.5R e (see Figure 1 of Saglia, Bender & Dressier 1993). Here (J) is the 



-4- 



average surface brightness inside R. If the galaxies considered are large (R e > 10 arcsec), no 
seeing corrections are needed (see Saglia et al. 1993). 

Until the use of CCD detectors, differential luminosity profiles of galaxies were obtained 
largely from photographic plates. The procedure required to calibrate the nonlinear response 
of the plates and to digitize them is very involved. As a consequence, it was possible to de- 
rive accurate luminosity profiles or two-dimensional photometry only for a small number 
of galaxies (see, for example, de Vaucouleurs & Capaccioli 1979). Using this sort of data, 
Thomsen & Frandsen (1983) derived R e and Mtot for a set of brightest elliptical galaxies 
in clusters at redshifts < 0.15. They fit a two-dimensional R 1 ^ law convolved with the 
appropriate point-spread function and briefly investigated the systematic effects of sampling 
(pixel size), signal-to-noise ratio, and shape of the profile on the derived photometric quan- 
tities. Lauberts & Valentijn (1989) digitized and calibrated the blue and red plates of the 
ESO Quick Schmidt survey to derive the photometric parameters of a large set of southern 
galaxies. Here the total magnitudes are not corrected for extrapolation to infinity, but are 
defined as the integrated magnitude at the faintest measured surface brightness (beyond the 
25 B mag arcsec -2 isophote) for which the luminosity profile is monotonically decreasing. In 
addition, the catalogue gives the parameters derived by fitting a "generalized de Vaucouleurs 
law" (I = I exp(—(r/a) N ); compare to Eq. [Uj) to the surface brightness profiles. 

The last 15 years have seen the increased use of CCDs for photometry. CCDs are 
linear over a large dynamic range, can be flatfielded to better than 1% and allow one to 
eliminate possible foreground objects during the analysis of the data. Large samples of CCD 
luminosity profiles for early-type galaxies have been collected by Djorgovski (1985), Lauer 
(1985), Bender, Ddbereiner & Mollenhoff (1988), Peletier et al. (1990), Lucey et al. (1991), 
J0rgensen, Franx & Kjaergaard (1995). Using CCDs one can derive photometric parameters 
by fitting a curve-of-growth to the integrated surface brightness profile. One major concern 
of CCD photometry is sky subtraction. If the CCD field is not large enough compared to 
the half-luminosity radii of the galaxies, then the sky value determined from the frame may 
be systematically overestimated (due to contamination of the sky regions by galaxy light), 
leading to systematically underestimated R e and Mtot- This problem might however be 
solved with the construction of very large chips or mosaics of CCDs (see MacGillivray et al. 
1993, Metzger, Luppino & Miyazaki 1995). 

Among the most recent studies of galaxy photometry is the Medium Deep Survey per- 
formed with the Hubble Space Telescope. Casertano et al. (1995) analyse 112 random 
fields observed with the Hubble Space Telescope Wide Field Camera prior to refurbishment 
to study the properties of faint galaxies. They construct an algorithm which fits the two- 
dimensional matrix of data points to perform a disk/bulge classification. The i? 1//4 and 
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exponential components are convolved with the point spread function (psf) of the HST and 
Monte Carlo simulations are performed to test the results. Disk-bulge decomposition is at- 
tempted only for a few cases (see Windhorst et al. 1994), because the data are in general 
limited by the relatively low signal-to-noise and by the spatial resolution. 

In order to derive total magnitudes, galaxy photometry involves extrapolation of curves 
of growth to infinity, and therefore relies on fits to the galaxy luminosity profiles. Recently, 
Caon, Capaccioli & D'Onofrio (1993, hereafter CCO) and D'Onofrio, Capaccioli & Caon 
(1994) focused on the use of the R lj/4 law to fit the photometry of ellipticals. CCO find a 
correlation with the galaxy size and argue that if an i? 1//n law (Sersic 1968, see Eq. [I~6| ) is 
used to fit the luminosity profiles, then smaller galaxies (log_R e [kpc] < 0.5) are best fitted 
with exponents 1 < n < 4, while larger ones (logi? e [kpc] > 0.5) have n > 4. Half-light radii 
and total magnitudes derived using these results may differ strongly from those using R 1 ^ 
extrapolations. Finally, Graham et al. (1996) find that the extended shallow luminosity 
profiles of BCG are best fit by R l l n profiles with n > 4. 

To summarize, the EFAR collaboration has collected photoelectric and CCD photometry 
for 736 galaxies (see Colless et al. 1993 and Paper III), 31% of which appear to be spirals 
or barred objects. The remaining 69% can be subdivided in cD-like (8%), pure E (12 %) 
and mixed E/SO (49%); the precise meaning of these classifications is explained in detail in 
§3.4 of Paper III. We derived circularly averaged luminosity profiles for all of the objects. 
Isophote shape analysis can only be reliably performed for the subset of our objects which 
are large and bright enough, and will be discussed in a future paper. Since 96% of the EFAR 
galaxies have ellipticities smaller than 0.4, the use of circularized profiles does not introduce 
systematic errors on the photometric parameters derived (see § |3.7| ) and has the advantage of 
giving robust results for even the smaller, fainter objects in our sample. The galaxies show a 
large variety of profile shapes. Typically, each object has been observed several times, using 
a range of telescopes, CCD detectors, and exposure times, under different atmospheric and 
seeing conditions, with different sky surface brightnesses. 

Deriving homogeneous photometric parameters from the large EFAR data set has re- 
quired the construction of a sophisticated algorithm to (i) optimally combine the multiple 
photoelectric and CCD data of each object, (ii) fit the resulting luminosity profiles with a 
model flexible enough to describe the observed variety of profiles, (iii) classify the galaxies 
morphologically and (iv) produce reliable magnitudes and half-luminosity radii. 

This paper describes our method as applied in Paper III. It explores the sources of 
random and systematic errors by means of Monte Carlo simulations, and develops a scheme 
to quantify the precision of the derived parameters objectively. The fitting algorithm searches 
for the best combination of the seeing-convolved, sky-corrected i? 1//4 and exponential laws. 
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This approach fulfils the requirement (ii) above: it produces convenient fits to the extended 
components of cD luminosity profiles, it models the profile range observed in E/SO galaxies 
(from galaxies with flat cores to clearly disk-dominated SOs), and it reproduces the surface 
brightness profiles of spirals. Moreover, for the E/SOs and spirals, this approach determines 
the parameters of their bulge and disk components, to assist classification (requirement (iii)). 
Finally, this approach minimizes extrapolation (requirement (iv)) which is the main source 
of uncertainty involved in the determination of magnitudes and half-luminosity radii. 

Would it be possible to reach the same goals with another choice of fitting functions? 
We demonstrate here (§3.8|) that the i? 1//n profiles quoted above can be seen as a "subset" of 
the i? 1 / 4 plus exponential models and therefore might not meet requirement (ii). In addition, 
for n > 4 they require large extrapolations and therefore might fail to meet requirement (iv). 
What is the physical interpretation of the two components of our fitting function? There are 
cases (the above cited cD galaxies and the galaxies with cores) where our two-component 
approach provides a good fitting function, but the "disk-bulge" decomposition is not physical. 
However, we argue that the systematic deviations from a simple i? 1 / 4 law observed in the 
luminosity profiles of our early-type galaxies are the signature of a disk. We will investigate 
this question further in a future paper, where the isophote shape analysis of the largest and 
brightest galaxies in the sample will be presented. Would it be worth improving the present 
scheme by, for example, allowing for a third component (a second i? 1//4 or exponential) to be 
fit? This could produce better fits to barred galaxies or to galaxies with cores and extended 
shallow profiles. However, it is not clear that the systematic errors related to extrapolation 
and sky subtraction could be reduced. Summarizing, the solution adopted here fulfils our 
requirements (i)-(iv). 

This paper is organized as follows. §|2] describes the three-step fitting technique. This 
involves the algorithm for the combination of multiple profiles of the same object (§ |2.1| ), 
our two-component fitting technique with the additional option of sky fitting (§ |2.2| ), and 
the objective quality assessment of the derived parameters (§|2.3|). §|3] presents the results of 
the Monte Carlo simulations performed to test the fitting procedure and assess the precision 
of the derived photometric parameters. We explore a large region of the parameter space 
(R e B,h, D/B, T, see § |2.2j for a definition of the parameters) and test the performance of 



the fitting algorithm (§ |3.1| ). In § |3.2| we investigate the systematic effects introduced by 
possible errors on sky subtraction and test the algorithm to correct for this effect (see §2.2|) . 
The influence of the limited radial extent of the profiles ( §3.3j ), of the signal-to-noise ratio 
(§ |3.4| ), and of seeing and pixellation (§ ^.5| ) are also investigated. The profile combination 
algorithm is tested in § |3.6| . In § |3.7| we assess the effectiveness of using the fitting algorithm 
to derive the parameters of bulge and disk components of a simulated galaxy. A number 



of different profiles are considered in § |3.8| to test their systematic effect on the photometric 
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parameters. We show that the R 1/>n profiles can be reproduced by a sequence of _R 1//4 plus 
exponential profiles, with small systematic differences (< 0.2 mag arcsec~ 2 ) over the radial 
range R e /20 < R < 5R e (see discussion above). In §[JH5] we discuss how to estimate the 
precision of the derived photometric parameters. In §[| we summarize our results in terms 
of the expected uncertainties on the derived photometric parameters. 

2. The fitting procedure 

The algorithm devised to fit the luminosity profiles of EFAR galaxies (see Paper III) in- 
volves three connected steps, (i) the combination of multiple profiles, (ii) the two-component 
fitting, and (iii) the quality estimate of the results. In the first step, the multiple CCD 
luminosity profiles available for each object are combined taking into account differences in 
sensitivity or exposure time, and sky subtraction errors. A set of multiplicative and addi- 
tive constants is determined (ki, Aj), which describe respectively the relative scaling due 
to sensitivity and exposure time and the relative difference in sky subtraction errors. The 
absolute value of the scaling is the absolute photometric calibration of the images. This is 
accomplished as described in Paper III, making use of the photoelectric aperture magnitudes 
and absolute CCD calibrations. The absolute value of the sky correction can be fixed either 
to zero or to a percentage of the mean sky, or passed to the second step to be determined 
as a result of the fitting scheme. 

The second step fits these combined profiles. The backbone of the fitting algorithm 
is the sum of the seeing- convolved i? 1//4 and the exponential laws. We have discussed the 
advantages of this choice in the Introduction. This combination produces a variety of lu- 
minosity profiles which can fit a large number of realistic profiles to high accuracy. The 
photometric parameters derived from this approach do not require large extrapolations, if 
the available profiles extend to at least 4R e . When galaxies with disk and bulge components 
(E/SOs and spirals) are seen at moderate inclination angles (as it is the case for the EFAR 
sample, where 96% of galaxies have ellipticities less than 0.4, see Paper III), then the algo- 
rithm is also able, to some extent, to determine the parameters of the two components. In 
Paper III this information is used, together with the visual inspection of the images and, 
sometimes, the spectroscopic data, to classify each EFAR object as E, E/SO, or spiral. While 
we believe that in these cases the two components of the fits are indicative of the presence 
of two physical components, additional investigation is certainly required to test this con- 
clusion. This will involve the isophote shapes analysis (Scorza & Bender 1995), the fitting 
of the two-dimensional photometry (Byun & Freeman 1995, de Jong 1996), the colors and 
metallicities (Bender & Paquet 1995), and the kinematics (Bender, Saglia & Gerhard 1994) 
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of the objects. We intend to address some of these issues in future papers for a selection of 
large and bright EFAR galaxies. 

The third step assigns quality parameters to the derived photometric parameters. Sev- 
eral factors determine how accurate these parameters can be expected to be. explores 
in detail the effects of sky subtraction errors, radial extent, signal-to-noise ratio, seeing and 
sampling, and goodness of fit. A global quality parameter based on these results quantifies 
the precision of the final results. 



2.1. Profile combination 

The first step of the fitting algorithm is to combine the multiple profiles available for 
each galaxy. Fitting each profile separately, and averaging the results produces severely 
biased results if the fitted profiles differ in their signal-to-noise ratio, seeing and sampling, 
radial extent, and sky subtraction errors. Only a simultaneous fit can minimize the biasing 



effects of these factors (see § |3.6[) . 



Apart from the very central regions of galaxies, where seeing and pixel size effects can 
be important, the profiles of the same object taken with different telescopes and instruments 
differ by a normalization (or multiplicative constant) only and an additive constant. The first 
takes into account differences in the efficiency and transparency, while the second adjusts for 
the relative errors in the sky subtraction. Let Ii(R), i — 1 to n denote the n available profiles 
in counts per arcsec 2 at a distance R from the center, and consider the profile I max {R) as 
the one having the maximum radial extent. In general the radial grids on which the profiles 
Ii(R) have been measured will not be the same, but it will always be possible to (spline) 
interpolate the values of I ma x{R) on each of the grid points of the other profiles I%{R). The 
normalization of the profiles relative to the profile I max (R) and the quantity A, 

(related to Aj/fcj the correction to the sky value of the profile Ii{R)) are the multiplicative 
and additive constants to be sought, so that: 

I[{R) = hh(R) - A,. (1) 

The ki and Aj constants are determined by minimizing the x 2 Tike functions (see the related 
discussion for Eq. |8|): 

x- = E w i( R ) (W^) - hh(R) + A t ) 2 • (2) 

R>R C 



The inner cutoff radius R c is 6 arcsec or half of the maximum extent of the profile, if this is 
less than 6 arcsec. This cutoff minimizes the influence of seeing, while retaining a reasonable 
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number of points in the sums. Here Wi(R) = l/<Ji(R) 2 are the relative weights of the data 
points, which are related to the expected errors for the profile If 



iGJiW + GiSkyi + RONf/S? 

ai{R) = , 3 

12-KR/Si 



where Si, Gi and RONi are the scale (in arcsec/pixel), the gain and the readout noise of 
the CCD used to obtain the profile /j (see Table 2 of Paper III). The denominator of Eq. 
[3] assumes that all of the pixels in the annulus at radius R ^ have been averaged to get 
I{R) and therefore underestimates the errors if some pixels have been masked to delete 
background or foreground objects superimposed on the program galaxies. If R = (i.e., the 
central pixel) the following equation is used: 

<n(R = 0) = sJg~UR = 0) + G.Sky, + RON? /Si (4) 

The weight in this fit monotonically increases with radius. The errors er Mi on the surface 
brightness magnitudes = — 2.5 log 1; are related to Eqs. and HI through: 



2.hai{R) loge 
h{R) 

By requiring d\\ / dki = and dxj / <9A j = we solve the linear system for ki and A j . 



At this stage the relative sky corrections are known for all of the profiles except the 
most extended one. This last correction A max can either be computed as part of the fitting 
program (see Eqs. |11] and [12]), or fixed to a given value. 



In §|3] the strategy of setting the mean percentage sky errors (for a given galaxy) to zero 
will be tested extensively against the above. For this case one requires: 

A <- + £_^ = . (6) 



Sky maa: ■ hSkYi 

In general, Eq. |6]is not a good choice and gives rise to systematic errors (see Fig. |6|), however 
it is preferred when the sky fitting solution (Eq. |12[) requires excessively large extrapolations. 
Forty percent of the fits presented in Paper III are performed using Eq. |6]. 

Note that for both Eq. || and described below, the value of A max is determined 
iteratively, by minimizing Eq. ||, having redefined I max (R) as I' max (R), where I' max (R) = 
Imax(R) — A max , and repeating the procedure until it convergences. Four or five iterations 
are needed to reach a precision < 10~ 5 when Eq. |] is used. Convergence is reached while 



performing the non-linear fitting of §[2.2|, when using Eq. O. Sky corrections, as computed 



in Paper III, are less than 1 % for 80% of the cases examined. 
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The absolute scaling, k max , of the I max (R) profile represents the photometric calibration 
of the profiles. This is performed as described in Paper III using the photoelectric aperture 
magnitudes and CCD zero-points. In the following we set k max = 1. 

2.2. i? 1 ' 4 -(-exponential law fitting 

The surface brightness profiles of each galaxy are modeled simultaneously by assuming 
that they can be represented by the sum of a de Vaucouleurs law (the "bulge" component 
indicated by B) and an exponential component (the "disk" component indicated by D): 

f(R, R eB , h, D/B, T, S) B+D = f B + f D , (7) 

where R eB is the half-luminosity radius of the bulge component, h the exponential scale 
length of the disk component, D/B the disk to bulge ratio, V the FWHM of the seeing profile, 
and S the pixel size. Both laws are seeing-convolved as described by Saglia et al. (1993) 
and take into account the effects of finite pixel size. Definitions and numerical details can be 
found in the Appendix. The results presented in Paper III show that Eq. [7] gives fits with 
respectably small residuals. The differences in surface brightness A/x = fi — fifa are typically 
less than 0.05 mag arcsec -2 , while those between the integrated aperture magnitudes are a 
factor two smaller. However our formal values of reduced \ 2 ( see discussion below) indicate 
that very few galaxies (less than 10%) have luminosity profiles that are fit well by the model 
disk and bulge. Over 90% of the fits have reduced x 2 larger than 2. In this sense Eq. |7] is 
not a statistically good representation of the galaxy profiles. 

A hybrid non-linear least squares algorithm is used to find the R eB , h, D/B and the 
vector of seeing values which gives the best representation f B+B) (R) of the profiles Ii(R), 
taking into account the sky corrections Aj/fcj. The algorithm uses the Levenberg-Marquardt 
search (Press et al. 1986), repeated several times starting from randomly scattered initial 
values of the parameters. The search is repeated using the Simplex algorithm (Press et al. 
1986). The best of the two solutions found is finally chosen. This approach minimizes the 
biasing influence of the possible presence of several nearly-equivalent minima of Eq. [|, a 
problem present especially when low disk-to-bulge ratios are considered (see discussion in 

All of the profiles Ii(R) available for a given galaxy are fitted simultaneously determining 
the appropriate value of the seeing Tj, for each single profile i. The minimization is performed 
on the function: 

xl b+d =j:( e n+ e n), (8) 

« \RXfB+D>-&i/ki R,Xif B +D<-Ai/ki ) 
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where: 



and: 



T> = -2.5 log 



T < = -2.5 lot 



XifB+pjR, ReB, h, D/B, Tj, Sj) + A % /k % 
k{R) 

KfB+D[R, ReB, h, D I B, Tj, Si 



7T> (9) 



Ii(R) — Ai/ki 



^- (10) 



The penalty function pi is introduced to avoid unphysical solutions and increases XtotB+D 
to very large values when D/B < or when the values of of R e B or h become too large 
(> 300") or too small (< 1"). The use of the T> and T< terms ensures that the arguments 
of the logarithm are always positive. The sky correction is usually applied to the fitting 
function (see Eq. |9|). However, data points where Xifs+D + Ai/ki < (this may happen 
when a negative sky correction Ai/ki is applied) are included using Eq. |TD|, which applies 
the sky correction to the data points. Note also that Eq. |^ is the weighted sum of the 
squared magnitude residuals. This is to be preferred to the weighted sum of the squared 
linear residuals, which is dominated by the data points of the central parts of the galaxies. 

The model normalization relative to the profile Ii(R), Aj, is determined by requiring 
dxlJdXi = 0, where: 

xl = $><(i2) ~ KfB+D(R) - ^Ihf ■ (ii) 

R 

Note that the ratios X max / \ can in principle differ from the constants ki, because of (residual) 
seeing effects (see, e.g., R c in Eq. |2|) and systematic differences between model and fitted 
profiles. In fact, the differences are smaller than 8% in 85% of the fits performed with more 
than one profile (see Paper III). When a bulge-only or a two-component model is used, 
the total magnitude of the fitted galaxy, in units of the I max (R) profile, is computed as 
M TO t = — 2.51og(L B + L D ), where L B = X max R\ B (see Eq. \Kl\ , with this normalization 
one has I e B = X max / (7 .22tt)) is the luminosity of the bulge and Lp = {D/B)Lb is the 
luminosity of the disk. When a disk-only model is used, then Mtot = — 2. 5 log L^, where 
Ld = X max h 2 (see Eq. |A2| , with this normalization one has Iq = X max / (2ti)) . Note again 
that the photometric calibration of these magnitudes Mtot to apparent magnitudes ttlt is 
performed in Paper III using photoelectric aperture magnitudes and CCD zero-points. 

The sky correction to the profile I max can be set to a given value (zero for no sky 
correction, using Eq. ^ for zero mean percentage sky correction). Alternately, a fitted sky 
correction A max can be determined by additionally requiring dx\ max /dA max = 0, where: 

X\ max = W max{R) {Imax{R) ~ X max f'B+D{R) ~ A max ) 2 . (12) 
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If the resulting A max produces \ ma xfB+D + A m(M . < at any R, Eq. [K] is used to compute 
the corresponding contribution to Eq. |8|. When using Eq. [L2], the constants ki and Aj 
are computed again using I' max {R) = I m ax(R) ~ A max (see the profile combination iterative 
algorithm in § j2.1| ). The Monte Carlo simulations of §^ show that Eq. [T2| gives an unbiased 
estimate of the sky corrections when the Jb+d is a good model of the fitted profiles. Eq. || 
is to be preferred when large extrapolations are obtained; 60% of the fits presented in Paper 
III are performed using Eq. |12"| . 

One might use the equivalent of Eq. [12] for the profiles h{R) to compute the corrections 
Aj directly from the fit, without having to go through Eq. |2|. This would automatically take 
into account the seeing differences of the profiles. However, tests show that this approach 
does not produce the correct relative sky corrections between the profiles, if the fitting 
function does not describe the fitted profiles well. Finally, one might try deriving Aj and Aj 
by minimizing Eq. [| for these two additional parameters. The adopted solution, however, 
speeds up the CPU intensive, non-linear minimum search, since Aj and Aj are computed 
analytically. 

The fit is repeated using a pure de Vaucouleurs law (D/B=0) and a pure exponential 
law (B/D=0). In analogy with Eq. two other Xtot are considered for these fits, XtotB 
and XtotD- A (conservative) 3a significance test (see discussion after Eq. |14|) is performed 
to decide whether the addition of the second component improves the fit significantly. The 
bulge-only fit is taken if: 



XtotB 1 <■ 3 



XtotB+D \ N- 



ifree-- (13) 



"B+D 



The disk-only fit is taken if: 



XtotD 1 <c 3 



XtotB+D \ N^+d 



(14) 



The number of degrees of freedom of the R 1 ^ 4 plus exponential law fit is N^ D = 
Ndata — N sky — 3 — 2N pro f, where N data is the number of data points involved in the sum of 
Eq. H N sky = 1 if the sky fitting is activated, zero otherwise, and 3 + 2N pro f are the number 
of parameters fitted (R e B,h,D/B, Mtot, N pro f seeing values and N pro f — 1 normalization 
constants Aj, where N pro f is the number of fitted profiles). If the errors are gaussian, 
XtotB+D follows a x 2 distribution of degrees of freedom. If the bulge plus disk model is 

a good representation of the data, then the XtotB+D ~ ^b+d * n the mean, with an expected 



dispersion J 2Ng+ e D . In this case Eqs. [H| and |14| are a 3a significance test on the conservative 
side, meaning that one-component models are preferred, if two-component models do not 
improve the fit by more than 3a. In fact, Paper III shows that only 10% of the fits are 
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statistically "good" (Xmb+d ~ n b+d)- The median reduced x 2 , X 2 = Xmb+dI n b+d, 
is ~ 6, indicating the existence of statistically significant systematic deviations from the 
simple two-component models of Eq. [7]. Fortunately, the tests performed in §|3| show that 
reliable photometric parameters can be obtained even in these cases. Note that fits based 
on the R l l n profiles do not give better results: Graham et al. (1996) obtain reduced x 2 ~ 10 
for their sample of brightest cluster galaxies. Eqs. and Eq. [14] as applied in Paper III 
select a bulge-only fit in 14 % of the cases, and a disk-only fit in less than 1%. In the 85% 
of the cases when both components are used, the median value of the significant test is 16a, 
with significance larger than 5cr in 90% of the cases. In the following sections and plots we 
shall indicate the reduced x 2 with x 2 ■ 

Total magnitudes of galaxies are extrapolated values. In order to quantify the effect of 
the extrapolation, we also derive the percentage contribution to Mtot = — 2.51og(Ls + Lp) 
due to the extrapolated light beyond the radius R ma x of the last data point. In 80% of 
the galaxies examined in Paper III this extrapolation is less than 10%. The half-luminosity 
radius R e (and the D n diameter, see Paper III) of the best-fitting function is computed using 
Eqs. and |A4| , so that seeing effects are taken into account. Finally, the contamination 
of the sky due to galaxy light is estimated by computing the mean surface brightness in the 
annulus with radii R™ ax and 2R™ ,ax 1 where R™ ax is the radius of the last data point of the 
profile i. Galaxy light contamination is less than 0.5 % of the sky in 80% of the cases studied 
in Paper III. 



Using the appropriate seeing-convolved tables (see § |2.2|) , the fitting algorithm can also 
be used to fit a \I/ = 12 model (see description in Saglia et al. 1993 and the appendix here) 
plus exponential, or a smoothed R 1 / 4 law plus exponential. These additional fitting models 
are useful to study the effects of the central concentration and radial extent of galaxies (see 



2.3. Quality parameters 

The third step in the fitting procedure assigns quality estimates to the derived photo- 
metric parameters. Several factors determine their expected accuracy, (i) Low signal-to-noise 
images provide fits with large random errors, (ii) Images of small galaxies observed under 
poor seeing conditions and/or with inadequate sampling (a detector with large pixel size) 
give systematically biased fits, (iii) Images of large galaxies taken with a small detector give 
profiles with too little radial extent and fits involving large, uncertain extrapolations, (iv) 
Sky subtraction errors bias the faint end of the luminosity profiles and therefore the fitted 
parameters. Finally, (v) bad fits to the luminosity profiles provide biased quantities. The 
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effects of these possible sources of errors are estimated by means of Monte Carlo simulations 
in §|. 

Based on these results, one can assign the quality estimates Q ma x, Qr, Qs/n, Qsky 
C^Sky' Q E) Q x 2 according to the rules listed in Table |], where increasing values of the 
quality estimates correspond to decreasing expected precision of the photometric parameters 
derived from the fits. The global quality parameter Q: 

Q = Max(Q max , Q r , Qs/n, Qsky ^Sky Qe, Q x 2 )> (15) 

assumes values 1,2,3, corresponding to expected precisions on total magnitudes AM TO t ~ 
0.05, 0.15, 0.4, on the logarithm of the half-luminosity radius A logi? e ~ 0.04, 0.1, 0.3 and on 
the combined quantity FP = logR e - 0.3(SB e ) AFP « 0.005, 0.01, 0.03 (see §g]| Fig. [TJ). 
Paper III shows that 16% of EFAR galaxies have Q = 1, 73% have Q = 2 and 11% Q = 3. 

Note that M TOT and, therefore, FP are subject to the additional uncertainty due to 
the photometric zero-point. In Paper III we extensively discuss this source of error and find 
that it is smaller than 0.03 mag per object, for all of the cases (86%) where a photoelectric 
or a CCD calibration has been collected. 



Table 1: The definition of the quality parameters 



Rmax / -F^e 


Qmax 


Rf/Tf 


Qr 


S/N 


Qs/N 


Extrap 


Qe 


< 1 

> 1,< 2 
> 2 


3 
2 
1 


< 2 
> 2 


2 
1 


< 300 
> 300 


2 
1 


> 0.3 
< 0.3 


3 
1 




x 2 


Q x 2 


^Sky - ( SB l ) 


^Sky 


|5Sky/Sky| 


^5Sky 






> 25 
> 12.5, < 25 
< 12.5 


3 
2 
1 


< 0.75 
> 0.75 


2 
1 


> 0.03 
> 0.01, < 0.03 
< 0.01 


3 
2 
1 







3. Monte Carlo Simulations 

The fitting procedure described in the previous section has been extensively tested on 
simulated profiles with the goals of checking the minimization algorithm and quantifying the 
effects of the errors described in §2.3j . Luminosity profiles of models with known parameters 
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have been fitted, to compare input and output values. In all of the following figures the 
output parameters of the fit are indicated with the superscript / for "fit" (for example, T^). 

As a first step (§|3j]-|3l| Figures PHl2l) , we ignore possible systematic differences between 
test profiles and fitting functions (such as the ones possibly present when fitting real galaxies, 
see discussion in Paper III) and generate a number of i? 1 / 4 plus exponential model profiles 
of specified R e B, h, D/B ratio, seeing F and total magnitude, using the seeing-convolved 
tables described in the Appendix. A constant can be added (subtracted) to simulate an 
underestimated (overestimated) sky subtraction. Given the pixel size, the sky value, the 
gain and readout noise, appropriate gaussian noise is added to the model profile following 
Eqs. |3] and The maximum extent of the profiles can be specified to simulate the finite 
size of the CCD. The profile is truncated at the radius where noise (or the sky subtraction 
error) generates negative counts for the model. The signal-to-noise ratios computed in the 
following refer to the total number of counts in the model profile out to this radius. The 
parameter space explored in all of the simulations discussed in § |3TT| - |3~5| is displayed in Figure 
H] and covers the region where the EFAR galaxies are expected to reside (see Paper III). 
Different symbols identify the models (see caption of Fig. |l|). As a second step (§ |3~T1 - |3.8| , 
Figures |13l-|l"5|), we explore the influence of systematic differences between test profiles and 
fitting functions. In § |3.7| we show that fitting circularized profiles of moderately flattened 
galaxies (as the ones observed in Paper III) allows good determinations of the photometric 
parameters and also of the bulge and disk components. In § |3.8j we fit the R l / n profiles, 
achieving two results. First, we quantify the influence of the quoted systematic effects on 
the fitted photometric parameters. Second, we suggest that the possible correlation between 
galaxy sizes and exponent n (see discussion in the Introduction) reflects the presence of a 
disk component in early-type galaxies. § |3.9| summarizes the results by calibrating the quality 
parameter Q of Eq. |15[ 



3.1. The parameter space 

In this section we discuss the results obtained by fitting the models indicated by the 
crosses in Figure [1]. For clarity the parameters are also given in Table |2|. No sky subtraction 
errors are introduced and the sky correction algorithm is not used. The detailed analysis 



of the possible sources of systematic errors discussed in § |3.2| - |3.4| is performed on the same 
sample of models. More extreme values of the parameters are used when testing the effects 
of seeing and resolution (§575) . The profiles tested in this section extend out to 4R e , have a 



pixel size of 0.4 arcsec and normalization of 10 7 counts, with Gi = RONi = 1 (see Eq. 
corresponding to S/N « 1000. 
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Table 2: The parameters of the models indicated by the crosses of Figure [l] (see § |3.1| ). A 
model for each combination of parameters in the two blocks separately has been generated. 
D/B = oo indicates exponential models {B/D = 0). 



Parameter 


Values 


ReB{ ) 


4, 8, 12, 16, 20, 32 


K) 


4, 8, 12, 16, 20, 32 


D/B 


0, 0.1, 1, 1.2, 1.6, 2, 3.2, 5, oo 


r(") 


1.5, 2.5 


Sky /pixel 


1000 


ReB{") 


2,3 


h(") 


3, 6 


D/B 


0, 0.2, 0.4, 0.8, 1.6, oo 


TO 


1.5, 2.5 


Sky/pixel 


500 
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Figure shows the precision of the reconstructed parameters. Total magnitudes are 
derived with a typical accuracy of 0.01 mag, R e and T to 3%, R £ b and h to ~ 8%, D/B to 
Rj 10%. The errors AM TO t = M TOT — Mf or and Ai? e = log R e /R{ are highly correlated, 
with insignificant differences from the relation AFP = AR e — 0.3(AMtot + 5AR e ) = 
AR e — 0.3A(SB e ). Galaxies with faint {D/B < 0.3) and shallow (h/R e B > 2) disks show the 
largest deviations. This partly reflects a residual (minimal) inability of the fitting program 
to converge to the real minimum x 2 (there are 3 points with x 2 > 10), but stems also from 
the degeneracy of the Bulge plus Disk fitting. Figure § (a) and (b) show the example of a 
model with D/B = 0.1 and h/R e B = 5 where a very good fit is obtained (x 2 = 1.3) yet 
there is a 0.05 mag error on Mtot and the disk solution is significantly different from the 
input model. Note that the largest deviations AMtot and AR e are associated with the 
largest extrapolations (~ 20%). In the following sections we shall see that extrapolation 
is the main source of uncertainty, when determining total magnitudes and half-luminosity 
radii. The uncertainties AR £ b on the bulge scale length are smallest with bright bulges, 
while those on the disk scale length Ah are smallest with bright disks. The algorithm to 
opt for one-component best-fits (Eqs. [13|-[14]) identifies successfully all of the one-component 
models tested (bulges plotted at \ogD/B = —1.1 and log h/R eB = —11, disks plotted at 
\ogD/B = 1.1 and \ogh/R e s = 1.1 in Figure |2|). For only two models (with D/B = 0.1 
and large h/R e B) is the bulge-only fit preferred (using the 3<r test) to the two-component fit 
(circled points in Figure |2|). 

Figure f| shows the results obtained by fitting a pure bulge or a pure disk. As before, 
no sky subtraction error is introduced and the sky correction algorithm is not activated. 
Neglecting one of the two components strongly biases the derived total magnitudes and 
half-luminosity radii. In the case of the R 1 ^ fits, already test models with values of D/B 
as small as ~ 0.2 give fitted magnitudes wrong by 0.2 mag, and R e by more than 30%. 
The systematic differences correlate with the amount of extrapolation involved, and large 
extrapolations yield strongly overestimated magnitudes and half-luminosity radii. However, 
the resulting correlated errors \AFP\ are almost always smaller than 0.03. In the case of pure 
exponential fits, the derived total magnitudes and half-luminosity radii are always smaller 
than the true values, since very little extrapolation (< 1%) is involved. Consequently, a 
positive, correlated error AFP (« +0.03) is obtained. Finally, note that pure bulge fits are 
bad fits of the surface brightness profiles (x 2 > 10), but may appear to give acceptable fits 
of the integrated magnitude profiles. (One can easily show that the differences in integrated 
magnitudes are the weighted mean of the differences in surface brightness magnitudes). 
Figure |3] (c) and (d) shows such an example for an R 1 / 4 fit to a model with D/B = 0.8 
and h/R e B = 1. The residuals in the integrated magnitude profile are always smaller than 
0.07 mag, but a x 2 = 181 is derived, with AM TOT = 0.32 and R{/R e = 1.65. These 
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considerations suggest that magnitudes and half-luminosity radii derived by fitting the R l / A 
curve of growth to integrated magnitude profiles (Burstein et al. 1987, Lucey et al. 1991, 
J0rgensen et al. 1995, Graham 1996) may be subject to systematic biases, as indeed Burstein 
et al. (1987) warn in their Appendix. This might be important for the sample of J0rgensen 
et al. (1995), where substantial disks are detected in a large fraction of the galaxies by means 
of an isophote shape analysis. It is certainly very important for the sample of cD galaxies 
studied by Graham (1996, see discussion in Paper III). These objects have luminosity profiles 
which differ strongly from an i? 1//4 law. 

Finally, note that the systematic errors shown in Figure f| (and in the figures of the 
following sections) cannot be simply estimated by considering the shape of the x 2 function 
near the minimum. Figure [5] shows the 1, 2, 3, 5cr contours of constant \ 2 f° r an R 1 ^ fit 
to a h/R e B = 0.5, D/B = 0.1 R l / A plus exponential model. The reduced x 2 (8.47 at 
the minimum) has been normalized to 1, so that la corresponds to a (normalized) \ 2 = 
1 + yj2/Nf ree = 1.11. The errors, estimated at the 5a contour, underestimate the differences 
between the fit and the model by a factor 2. This results from the extrapolation involved 
and can be as large as one order of magnitude for models with larger D/B ratios. 

3.2. Sky subtraction errors 

Sky subtraction errors can induce severe systematic errors on the derived photometric 
parameters of galaxies. Figure § shows the parameters derived from the i? 1 / 4 plus exponential 
models examined in the previous section, where now the sky has been overestimated or 
underestimated by ±1%. The sky correction algorithm is not activated. 

The biases become increasingly large as the sky brightness approaches the effective 
surface brightness of the models. As expected, underestimating the sky (a negative sky 
error) produces total magnitudes that are too bright and half-luminosity radii that are too 
large relative to the true ones. The size of the bias correlates with the extrapolation needed to 
derive M^ OT . The opposite happens when the sky is overestimated, but the amplitude of the 
bias is smaller, because there is no extrapolation. The correlated error AFP remains small 
(~ 0.05), except for the cases where large extrapolations are involved. The D/B ratio is ill 
determined, with better precision for models with extended disks (h/R e B > 2.5). The scale 
length of the bulge is better determined for low values of D/B (dominant bulge), the scale 
length of the disk component is better determined for large values of D/B (dominant disk). 
The parameter least affected is the value T of the seeing, which is determined in the inner, 
bright parts of the models, where sky subtraction errors are unimportant. Bulge-only or disk- 
only models appear to be fit best by two-component models (crosses and triangular crosses 
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Fig. 1. — The parameter space of the _R 1//4 plus exponential profile of the Monte Carlo 
simulations discussed in Figure |2|-|l"l"|. Models of Figs. |2|-|7]: crosses (see also Table 2). Models 
of Fig. ||: skeletal triangles. Models of Fig. |9[ open triangles. Models of Figs. p7D| - p7T| : open 
squares. Models of Fig. [L2|: open pentagons. Models of Fig. [L3|: open hexagons. The small 
dots show the position of the EFAR galaxies as determined in Paper 3. The parameters of 
bulge only models are shown with h = 0. The parameters of disk only models are shown at 
R e B = 0. See discussion in 
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Fig. 2. — The reconstructed parameter space for the models indicated by the crosses in Fig. 
|]. No sky error is present. The quantities plotted on the y-axis are defined as AMtot — 
M tot -M£ ot , AR e = \ogR e /R f e , AFP = AR e -0.3(AM TOT + 5AR e ) = AR e -Q.3A(SB e ), 
A(D/B) = \og[(D/B)/(D/B)f), AR eB = \ogR eB /R f eB , Ah = logh/h?, AT = logiyr'. On 
the x-axis, the first three boxes show the input parameters of the models in the logarithm 
units (\ogD/B, log h/R e B, log R e /T). The last three boxes show the differences in mag- 
nitudes between the assumed sky value and the average effective surface brightness of the 
models (jWgky ~~ (SBe)), the logarithm of the reduced x 2 ; an d the fraction of light extrapo- 
lated beyond R ma x used in the determination of M^ OT . Models with D/B = (pure i? 1//4 
laws) are plotted at \ogD/B — — 1.1 and log h/R eB = —1.1. Models with B/D = (pure 
exponential laws) are plotted at logD/B = 1.1 and log/i/ ' R eB = 1.1. Models with D/B ^ 



which have been fitted with one component are circled. See §3.1| for a discussion of the 
results. 
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Fig. 3. — (a) A circular disk plus bulge model with D/B = 0.1 and h/R e B = 5 (crosses). 
The dotted curves show the luminosity profiles fi(R) = — 2.5 log I (R) of the bulge and the 
disk components, the dashed curve the fitted disk component, (b) The differences A/x (in 
mag arcsec -2 , open squares) between model surface brightness and the fitted one (dotted 
curve) (see Q. (c) The R l ' A fit (solid curve) to the surface brightness magnitude profile of 
a circular disk plus bulge model with D / B = 0.8 and h/R e B = 1 (crosses, one point in every 
four), (d) The differences Amag between the R 1 ^ integrated magnitudes and the fitted ones 
(solid curve, see §3.1| ). Note that \Amag\ < 0.07 even if large deviations A/x are present. 



(e) The fit to the circularized profile of a flattened bulge plus an inclined disk model (see 
§. |3.7| ). The luminosity profile of the model (crosses, one point in every four; the bulge and 
the disk components, with the listed parameters, are the full curves) is best fitted by an 
R 1 / 4 plus exponential law (dashed curves) with parameters R[ e = 16.34 arcsec, = 13.93 
arcsec, (D/B)f = 0.13, R e = 17.51 arcsec. (f) The residuals A/x of the fit (open squares, 
one point in every four) and the differences between the growth curves Amag (full curve). 
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Fig. 4. — The effects of fitting disk plus bulge test profiles by either a single bulge (first 
three rows of plots) or a single disk (last three rows of plots) model. The test models are 
indicated by the crosses of Fig. 0. AMtot, Ai? e , and AFP are defined as in Fig. 0, x-axis 
as in Fig. 0. See § p.l| for the a discussion of the results. Note the change of scale on the 
ordinate axis with respect to Figure |2|. 
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Fig. 5. — Illustration of the underestimation of the errors. The contours of constant \ 2 near 
the minimum of an R 1 ^ fit to a h/R e B = 0.5, D/B = 0.1 disk plus bulge model. The cross 
shows the best-fit solution, the circle near the upper left corner gives the real parameters of 
the model. The errors estimated at the 5cr contour underestimate the differences between 



the model and the fit by a factor 2 (see § |3~T 
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in Figure ^j. Finally, note that reasonably good fits (x 2 < 10) to the surface brightness 
profiles are always obtained, in spite of the large errors on the reconstructed parameters. 

The biases discussed above can be fully corrected when the sky-fitting algorithm of Eq. 
12| is applied. Figure [7| shows the reconstructed parameters of the models considered in § |3.1| , 



where sky subtraction errors of 0, ±1%, ±3% have been introduced. For most of the models 
examined, the errors on the derived quantities are no more than a factor 2 larger than those 
shown in Fig. 0. The sky corrections are computed to better than 0.5% precision. Larger 
errors AMtot and AR e are obtained for models with relatively weak {D/B < 0.3) and 
extended disks (h/R e B > 2.5), where the degeneracy discussed in §[D] is complicated by the 
sky subtraction correction. These cases give reasonably good fits (x 2 < 10), but are identified 
by the large extrapolation (> 0.3) involved. Models with concentrated disks {h/R eB < 0.2) 
can also be difficult to reconstruct, when h/Y ~ 1. For some of these problematic fits, 
one-component solutions are preferred by Eqs. |13| - |T^ (circles in Eq. |^). 

A common problem of CCD galaxy photometry is the relatively small field of view, 
particularly with the older smaller CCDs. If the size (projected on the sky) of the CCD is 
not large enough compared to the half-luminosity radius of the imaged galaxy, then the sky 
as determined on the same frame will be contaminated by galaxy light and biased to values 
larger than the true one. Total magnitudes and half-luminosity radii can therefore be biased 
to smaller values, the effect being more important for intrinsically large galaxies, which tend 
to have low effective surface brightnesses. The mean surface brightness in the annulus with 
radii R™ ax and 2R™ ax (see § |2.2|) predicted by the fit allows us to estimate the size of the 
contamination. 



3.3. Radial extent 

Photoelectric photometry of large, nearby galaxies rarely goes beyond 1 or 2 R e (Burstein 
et al. 1987) and the same applies for the surface photometry obtained with smallish CCDs. 
The typical profiles obtained in Paper III extend to a least 4 R e , but a small fraction of them 
are less deep, reaching 1 or 2 R e only. Here we investigate the effect of the radial extent of 
the profiles, keeping the normalization of the profiles fixed (10 7 counts, S/N ~ 10 3 ). Sky 
subtraction errors of 0, ±3% are introduced and the sky fitting is activated. Figure ^| shows 
the cumulative distributions of the errors on the derived photometric parameters as derived 
from the simulations, for a range of R max values. When R max = Re, rather large errors are 
possible (0.3 mag in the total mag nitude, > 30% in R f e ). The main source of error is again 
the large extrapolation involved when R max ~ R e , coupled with the sky correction which 
becomes unreliable for these short radial extents. As soon as R ma x > 3i? e the errors reduce 
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Fig. 6. — The biases introduced by a ±1% sky subtraction error. Quantities plotted as in 
Fig. 0. Models with D/B = which have been fitted with two components are shown as 
crosses. Models with B/D = which have been fitted with two components are shown as 
triangular crosses. Note the change of scale on the ordinate axis with respect to Figure 0. 
See discussion in 5R72". 
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Fig. 7. — The effects of the sky fitting algorithm. The parameters of the models of Fig. 
|6] with the sky subtraction errors of 0, ±1%, ±3%, are reconstructed using the sky fitting 
algorithm. Quantities and symbols plotted as in Figures || and |j. In addition, the difference 
AdSky=dSky/Sky-dSky^/Sky on the sky correction is plotted. Note the change of scale on 
the ordinate axis with respect to Figure ||. See discussion in § [3~2| . 
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to the ones discussed in §3.1| . The same kind of trend is observed for the parameters of 
the two component (A(D/B), AR e s, Ah). The seeing values are less affected, as they are 
sensitive to the central parts of the profiles only. Finally, note that in all cases very good 
fits are obtained (x 2 ~ 1). 



3.4. Signal-to-noise ratio 

For most of the galaxies discussed in Paper III, multiple profiles are available with 
integrated signal-to- noise ratios S/N > 300, the normalization used in the previous sections. 
But for some of the luminosity profiles a smaller number of total counts has been collected 
(see Figure |l[). Here we investigate how the signal-to- noise ratio of the profiles affects the 



outcome of the fits. As before, the subset of models of §^2] is used with R ma x < 4i? e (see 
comment at the beginning of §|3|). Sky subtraction errors of 0, ±3% are introduced and the 
sky fitting is activated. Figure ^ shows how the errors on the derived parameters increase 
when the signal-to- noise ratio is reduced. For fluxes as low as about 10 5 {S/N m 10 2 ) all of 
the derived photometric parameters become uncertain (0.2 mag in the total magnitudes, 20% 
variations in the derived R e , large spread A(D/B), AR eB , Ah, AT), as large extrapolations 
and uncertain sky corrections are applied. In all cases very good fits are obtained (x 2 ~ 1). 



3.5. Seeing and sampling effects 

Some of the galaxies considered in Paper III are rather small, with R e < 4". Here 
we investigate the effects of seeing and pixel sampling, when R e m F ^pixel size. Figure 



TU| shows that reliable parameters can be derived down to R e ~ T, with pixel sizes 0.4-0.8 



arcsec, with only a small increase of the scatter for R e < 2T. 

A small systematic effect is caused by the choice of the psf. Saglia et al. (1993) 
demonstrate that a good approximation of the psfs observed during the runs described in 
Paper III is given by the 7 psf with 7 = 1.5 — 1.7. We adopt 7 = 1.6 for the fits. Here 
we test the effect of having 7 = 1.5 or 1.7 with a pixel size of 0.8 arcsec. Figure [TTJ shows 
that if 7 = 1.5 is the true psf of the observations, then the half- luminosity radius, the total 
luminosity, the scale length of the bulge will be slightly overestimated, and the disk to bulge 
ratio slightly underestimated. A small systematic trend is observed in the correlated errors 
AFP. The scale length of the disk component is less affected. The sky corrections are also 
biased, but do not strongly affect the photometric parameters, because of the high average 
surface brightness of the small R e models. Seeing values suffer a very small, but systematic 
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Fig. 8. — The effect of the radial extent of the profiles on the precision of the derived pa- 
rameters. The cumulative distributions of the errors on the derived photometric parameters 
as derived from the simulations are shown for a range of R ma x values (full lines: R ma x — Re, 
dotted lines: R max = 2i? e , dashed lines R ma x = 3i? e , long-dashed lines: R ma x = 4i? e ). Good 
reconstructions are obtained when R max /R e > 2 (see §|3.3|). 
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Fig. 9. — The effect of the signal-to- noise ratio of the profiles on the precision of the derived 
parameters. Good reconstructions are obtained when S/N > 300 (see §|3~H|). Note the change 
of scale on the ordinate axis with respect to Figure ^| 
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effect. The opposite trends are observed if the true 7 is 1.7. In all cases very good fits are 
obtained (x 2 ~ 1). The systematic differences become unimportant for R e > 2T. 

Finally, the seeing values derived can be systematically biased, if the central concen- 
tration of the fitted galaxies does not match the one of the R 1 ^ plus exponential models. 
We investigate this effect by fitting the \l/ = 12 plus exponential or the smoothed i? 1 / 4 plus 



exponential models discussed in § [2.2| . We find that in the first case the seeing value is un- 
derestimated which compensates for the higher concentration of the \1/ = 12 component. 
The shallow radial decline of the luminosity profile in the outer parts introduces systematic 
biases in the reconstructed parameters, similar to those discussed for the R l ' n profiles, for 
large values of n (see §|3.8Q. The half-luminosity radii and total magnitudes derived are un- 
derestimated by 20% and 0.2 mag respectively, when a \l/ = 12 model with no exponential 
component is fitted. The biases are reduced when models with an exponential component 
are constructed. In the case of the smoothed i? 1//4 law, the seeing value is overestimated to 
fit the lower concentration of the smoothed i? 1 / 4 component. No biases are introduced on 
the other reconstructed parameters. 



3.6. Tests of profile combination 

In order to test the combination algorithm described in |2.1| , four profiles with different 
T, pixel sizes, normalizations, gain, readout noise, and sky subtraction errors (see Table || 
these parameters match the typical values of the profiles of Paper III) are generated for 
the set of models identified by the open pentagons of Figure [I]. Figure [12] shows the result 
of the test. The abscissa plots the residuals A of the parameters derived using the fitting 
procedure with profile combination. A dSky and Ar are averaged over the four obtained 
values. The ordinate plots the mean of the residuals of the parameters derived by fitting each 
single independently as crosses, and the residuals of each fit as dots. The profile combination 
algorithm obtains better precision on all of the parameters with the exception of T, where 
the maximum deviation is in any case smaller than 8%. 
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Fig. 10. — The effect of seeing and pixel sampling of the profiles on the precision of the 
derived parameters. Different symbols indicate different pixel sizes (small dot 0.4 arcsec, 
triangles 0.6 arcsec, squares 0.8 arcsec). Note the expanded ordinate scale with respect to 
Figures f§-|9. See discussion in §|5T5 . 
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Fig. 11. — The effect of the choice of the psf on the precision of the derived parameters. 
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Fig. 12. — The profile combination algorithm and the precision of the derived parameters. 
The x-axis plots the residuals A of the parameters derived using the fitting procedure with 
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plots the mean of the residuals of the parameters derived by fitting each single independently 
as crosses, and the residuals of each fit as dots (see discussion in §|3.6|). 
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Table 3: The parameters of the multiple profiles test (see §3.6|). 



Profile 


1 


2 


3 


4 


Pixel size (") 


0.4 


0.606 


0.862 


0.792 


Sky per pixel 


300 


350 


250 


1500 


<5Sky/Sky 


+1% 


-0.5% 


+ 1.5% 


+0.5% 


Rraax / Re 


4 


3 


4.5 


2.5 


Normalization 


10 7 


5 x 10 6 


10 7 


5 x 10 6 


Gain 


1 


3 


1 


2 


Ron 


1 


4 


1 


5 


r(") 


2 


2.1 


1.5 


2.4 



3.7. "Bulge" and "Disk" components 

The discussion of the previous sections shows that for a large fraction of the parameter 
space, i.e. when deep enough profiles are available, with large enough objects, not only 
can the global photometric parameters R e and Mtot be reconstructed with high accuracy, 
but also the parameters of the _R 1//4 and the exponential components. Here we investigate 
further if reliable "bulge" and "disk" parameters can be derived, when the profiles analysed 
are constructed from the superposition of these two components. 

With this purpose, we constructed a number of two-dimensional frames (filled triangles 
in Figure [I]) as the sum of a flattened i? 1//4 bulge and an exponential disk of given inclination. 
The bulge (disk) frames follow an exact i? 1//4 (exponential) law with R e B = Yl^jbja arcsec 

(h = 10^/cos(i) arcsec) along the minor axis. Three flattenings of the bulge (6/ a = 1, 0.7, 0.4), 
four inclinations for the disk (i = 0°, 30°, 60°, 80°, where i = 0° is face-on and i = 90° edge-on) 
and five values of the disk to bulge ratio (D/B — 0, 0.5, 1, 2, oo) are considered. The resulting 
models are normalized to 10 7 counts. The pixel size is 0.6 arcsec. The circularly averaged 
luminosity profiles are derived following the same procedure adopted for the observed galaxies 
(see Paper III) and extend out to ~ 4 — 6R e . A 1% sky error is introduced and the sky 
fitting procedure is activated. Note that the maximum flattening of the EFAR galaxies is 
b/a = 0.5, with 96% of the galaxies having b/a > 0.6 (see Paper III). This corresponds to 
(pure) disk inclinations i < 60°. 



Figure [1^ shows the reconstructed parameters as a function of the inclination angle 
of the disk, for the different flattenings of the bulge, using the sky fitting procedure. The 
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horizontal bars show models with D/ B = 0.5. The plot at the bottom right shows the scale 
lengths of the flattened bulge (filled symbols) or of the inclined disk as a function of the 
flattening angle (open symbols, % = arccos(b/a)) or of the inclination angle, normalized to 
the b/a = 1 or % = 0° values. When D/B is low (< 0.5) the errors are very small for every 
inclination angle. For larger values of D/B, reliable photometric parameters are obtained 
for i < 60°, but as soon as the disk is nearly edge-on, total magnitudes and half-luminosity 
radii are overestimated (by 0.1 mag and 20% respectively). The integrated circularized 
profiles, in fact, converge more slowly than the ones following the isophotes. The correlated 
errors AFP always remain very small. Similarly, the parameters of the two components 
are reconstructed well for i < 60°, but badly underestimate the disk when it is nearly edge- 
on. However, a decent fit is obtained, by increasing the half-luminosity radius of the bulge 
component (see Fig. ^ (e) and (f)). The sky correction is returned to better than 0.5% for 
i < 80°. The systematic effects connected to the flattening of the bulge are small for the 
range of ellipticities considered here (b/a > 0.6). 

These results indicate two potential problems, (i) galaxies may be misclassified due to 
the presence of an edge-on disk component not being recognized, or (ii) the photometric 
parameters may be systematically overestimated. However, these problems do not apply to 
the EFAR sample, where b/a > 0.5 always and b/a > 0.6 for 96% of the galaxies. Therefore, 
galaxies with bright edge-on disks are only a very small fraction. Galaxies with faint edge-on 
disks, which may not show large averaged flattenings, have low D/B ratios and therefore are 
not affected by problem (ii). In a future paper we will address the question whether in these 
cases the isophote shape analysis might detect these faints disks and improve on point (i). 

Finally, the two-dimensional frames described here have been used to calibrate the 



estimator of the galaxy light contamination described in §2.2j . We measured the sky in the 
same way as for the real frames of Paper III, by considering some small areas around the 
simulated galaxies. We find that the predicted galaxy light contamination overestimates the 
measured sky excess by at least a factor two, and therefore can be used as a rather robust 
upper limit to the galaxy light contamination. 



3.8. R l l n luminosity profiles 

The tests described above show that our fitting algorithm is able to reconstruct the 
parameters of a sum of an i? 1//4 plus an exponential law accurately. In these cases sky 
subtraction errors can also be corrected efficiently Even so, we do find in Paper III that 
luminosity profiles of real early-type galaxies show systematic differences from i? 1//4 plus 
exponential profiles, yielding to a median reduced x 2 °f 6- Here we quantify the systematic 
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Fig. 13. — The reconstructed parameters of the bulge plus disk models as a function of the 
inclination i of the disk. Different symbols indicate different flattenings of the bulge. The 
horizontal bars show models with D/B = 0.5. The plot at the bottom right shows the scale 
lengths of the flattened bulge (open symbols) or of the inclined disk (filled symbols) as a 
function of the flattening angle (i = arccos(b/a)) or of the inclination angle, normalized to 
the b/a = 1 or i = 0° values. Good reconstructions of the parameters are obtained when the 
inclination is less than 60° (see §|3.7|). 
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effects that would be produced in this case, by studying the case of the R x ' n profiles. 

CCO fitted the luminosity profiles of 52 early-type galaxies using the R 1 ^ law introduced 
by Sersic (1968): 

"(*r 



I(R) = 1210 



-1 



(16) 



where b n ~ 0.868n — 0.142, R™ is the half-luminosity radius, and 7™ the surface brightness 
at R™. The total luminosity is L T = KJ^Rf, where \ogK n « 0.03[log(n)] 2 + 0.441 log(n) + 
1.079. Eq. ^ reduces to Eq. [Al] for n = 4 and to Eq. |A2| for n = 1. For large values of n, 
Eq. [T^ describes a luminosity profile which is very peaked near the center and has a very 
shallow decline in the outer parts. Ciotti (1991) computes the curve of growth related to Eq. 
|16| analytically for integer values of n and finds that while already ~ 13% of the total light 
is included inside R < 0.05i2", only 80% of the total light is included inside 6R™ for n = 10. 

We fitted Eq. [161 modified to have a core at R < 0.05R2, to an i? 1 / 4 plus exponential 



model for n = 0.5 to n — 15 out to 6R™ . Fig. O shows the results of the fit for a selection of 



models. With the exception of the n = 0.5 model, all of the R x l n profiles can be described 
by a combination of an i? 1//4 and an exponential component, with residuals less than 0.2 mag 
arcsec -2 for R < 4R e . For n < 4 the residuals increase to 0.4 mag arcsec -2 at R > 5R™, 
where the fits are increasingly brighter than the R}l n profiles. For large values of n the 
residuals reach -0.4 mag arcsec -2 at R > 5-R™, where the fits are increasingly fainter than 
the R x l n profiles. The relation between n and the parameters of the decomposition is shown 
in Fig. pl| Models with 1 < n < 4 are fitted using a decreasing amount of the exponential 
component, with a scale length comparable to the one of the i? 1 / 4 component. Models with 
n > 4 are fitted with an increasing amount of the exponential component, with increasingly 
large scale length. Half-luminosity radii are progressively underestimated, being w 60% of 
the true values at n — 15. Correspondingly, total magnitudes are also underestimated, by 
0.25 magnitudes at n = 15. 

A possible problem can emerge for large values of n, if the sky fitting algorithm is 



activated. The dotted curves in Fig. 15 show that if the sky subtraction algorithm is 
activated (Eq. [12D, then larger systematic effects are produced. Note that the computed sky 
correction (dotted curve of Fig. |15|) is ~ for n w 1 or n w 4 only. For n > 4 the correction 
is used to reduce the systematic negative differences in the outer parts of the profiles. A 
comparison between the fitted sky corrections and the upper limits on the possible galaxy 



light contamination (see § |2.2| and |3.7|) gives an important consistency check. In the case 



shown in Fig. [15] the fitted sky corrections are twice as large as the upper limits on the 
galaxy light contamination. In a real case this, together with the rather large values of x 2 , 
would hint at an uncertain fitted sky correction. 
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The fact that the R 1 ^ sequence can be approximated by a subsample of R 1 ^ plus 
exponential models suggests a possible reinterpretation of CCO's results: the variety of 
profile shapes of early- type galaxies is caused by the presence of a disk component. Moreover, 
the use of the R}l n profiles to determine the photometric parameters of galaxies of large n is 
dangerous, since the extrapolation involved is large and the fitted profiles barely reach 2 or 
3-R™ , as derived from the fit. This problem is much smaller using the i? 1 / 4 plus exponential 
approach. 



3.9. Discussion 

We conclude our tests by discussing the quality parameters defined in Table [l] and their 
use to estimate the size of the systematic errors present. 

The definitions given in Table [I] have been derived after inspection of Figures |2| to 15, 
with the desired goal of identifying three classes of precision, AMtot < 0.05, AMtot < 0.15, 
AMtot > 0.15. The parameters Qe, Qmax, Q x 2 , Qs/n, and Qr are directly related to the 
simulations. Their low values imply that the fits involve a small extrapolation, extend to 
large enough radii, give low surface brightness residuals with a large enough signal-to-noise 
ratio and good spatial resolution. The definitions of Qg^y and Q^gky deal with the accuracy 
of the sky subtraction, taking into account that high surface brightness galaxies suffer less 
from this problem, and that large sky corrections indicate a lower quality of the data. Low 
values of Q (see Eq. |15|) imply low values of all quality parameters. 

Figure |TB| shows the cumulative distributions of the errors AMtot, Ai? e and AFP 
derived from all the performed disk plus bulge fits with sky correction algorithm activated, as 
a function of the different quality parameters. The two most important parameters regulating 
the precision of the photometric parameters are the level of extrapolation and the goodness 
of the fit, followed by the sky subtraction errors. A low Qe fixes the maximum possible 
overestimate of the parameters. A low Q x 2 with a low Qe constrains the underestimate 
and the reliability of the sky correction. The ranges of the errors match the desired goal of 
identify three classes of precisions. 

Finally, it is sobering to note that the constraints needed to achieve Q — 1, high 
precision total magnitudes and R e are rather stringent. Only 16% of EFAR galaxies have 
Q = 1. Most of the existing published values of Mtot and R e of galaxies are far below this 
precision, because of the restricted radial range probed by photoelectric measurements or 
small CCD chips, because of sky subtraction errors, and also by the use of the pure 
curve of growth fitting (see Figure f|). The related observational problems can be somewhat 
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Fig. 14. — The fits to the R lj/n law. Two plots are drawn for each value of n (given in the 
top right corner). In the top plot the crosses (one point in every seven) show the luminosity 
profiles /j,(R) = -2.5 log I(R) of the R 1/n law as a function of R/R e . The dotted and dashed 
curves show the best-fitting R 1 ^ and exponential laws respectively. In the bottom plot the 
residuals (full curves) in mag arcsec -2 from the fits to the R l / n law are shown. The dashed 
curve shows the residuals (in mag) from the curves of growth (see discussion in §[3] 
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Fig. 15. — The relation between n and the parameters of the decomposition (see § |3.8| ). The 
full curves refer to the results obtained with no sky subtraction errors. The dotted curves 
show the results obtained when the sky fitting algorithm is activated. 
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reduced with the use of large CCDs (see Introduction), but the a priori limiting factor of 
galaxy photometry, the extrapolation, will always remain with us at a certain level. 

On the other hand, the errors on M T ot and R e are strongly correlated, so that the 
quantity logi? e — 0.3{SB e ) is always well determined. This fact allows the accurate distance 
determinations achieved using the Fundamental Plane correlations despite the systematic 
errors in the photometric quantities. 

4. Conclusions 

We constructed an algorithm to fit the circularized profiles of the (early-type) galaxies 
of the EFAR project, using a sum of a seeing-convolved R 1 ^ and an exponential law. This 
choice allows us to fit the large variety of profiles exhibited by the EFAR galaxies homo- 
geneously. The procedure provides for an optimal combination of multiple profiles. A sky 
fitting option has been developed. A conservative upper limit to the sky contamination due 
to the light of the outer parts of the galaxies is estimated. From the tests described in 
previous sections we draw the following conclusions: 

1) The reconstruction algorithm applied to simulated R 1 ^ 4 plus exponential profiles 
shows that random errors are negligible if the total signal-to-noise ratio of the profiles exceeds 
300. Systematic errors due to the radial extent of the profiles are minimal if R max jR e > 2. 
Systematic errors due to sky subtraction are significant (easily larger than 0.2 mag in the 
total magnitude) when the sky surface brightness is of the order of the average effective 
surface brightness of the galaxy. They can be reliably corrected for as long as the fitted 
profiles show small systematic deviations (x 2 < 12.5). 

2) Strong systematic biases (errors larger than 0.2 mag in the total magnitudes) are 
present when a simple i? 1 / 4 or exponential model is used to fit test profiles with disk to 
bulge ratios as low as 0.2. 

3) The use of the shape of the (normalized) x 2 function badly underestimates the (sys- 
tematic) errors on the photometric parameters. 

4) Systematic biases emerge when test profiles are derived for systems with significant 
disk components seen nearly edge-on, or when the fitted luminosity profile declines more 
slowly than an i? 1 / 4 law. The parameters of bulge plus disk systems can be determined to 
better than « 20% if the disk is not very inclined (i < 60°). 

5) The sequence of R l ^ n profiles, recently used to fit the profiles of elliptical galaxies 
by Caon et al. (1993), is equivalent to a subset of R 1 ^ and exponential profiles, with 
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Fig. 16. — The precision of the reconstructed total magnitudes Mtot, the half-luminosity 
radii R e and the combined quantity FP = logi? e — 0.3(SB e ). The cumulative distributions 
of the errors AMtot, AR e and AFP derived from all the performed disk plus bulge fits with 
sky correction algorithm activated are shown as a function of the different quality parameters 
defined in §2.4 The full lines plot the distributions when the parameters have value of 1, 
the dotted ones when the value is 2, the dashed ones when the value is 3. The distributions 
derived by selecting on the global quality parameter Q match the precision ranges identified 
in 50. 
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appropriate scale lengths and disk-to-bulge ratios, with moderate systematic biases for n < 8 
and residuals less than 0.2 mag arcsec -2 for R < 4R e . This suggests that the variety of 
luminosity profiles shown by early-type galaxies is due to the frequent presence of a weak 
disk component. 

6) A set of quality parameters has been defined to control the precision of the estimated 
photometric parameters. They take into account the amount of extrapolation involved to 
derive the total magnitudes, the size of the sky correction, the average surface brightness of 
the galaxy relative to the sky, the radial extent of the profile, its signal-to-noise ratio, the 
seeing value and the reduced \ 2 of the fit. These are combined into a single quality parameter 
Q which correlates with the expected precision of the fits. Errors in total magnitudes Mtot 
less than 0.05 mag and in half-luminosity radii R e less than 10% are expected if Q = 1, and 
less than 0.15 mag and 25% if Q = 2. 

89% of the EFAR galaxies have fits with Q = 1 or Q = 2. The errors on the combined 
Fundamental Plane quantity FP = logi? e — 0.3 (SB e ), where (SB e ) is the average effective 
surface brightness, are smaller than 0.03 even if Q = 3. Thus systematic errors on Mtot 
and R e only marginally affect the distance estimates which involve FP. 
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A. The fitting function 

The fitting procedure described in §Q assumes that fitted profiles can be well represented 
by the sum of a de Vaucouleurs (1948) law of half-luminosity radius R £ b and a surface 
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brightness I eB at R eB (with B for bulge component): 

I B (R) = I eB exp{-7.67[(R/R eB y/ 4 - 1]}, (Al) 

and an exponential component with exponential scale-length h and central surface brightness 
I (with D for disk component): 

I D (R) = I exp(-R/h). (A2) 
The R l/A law curve-of-growth is: 



F B {R) 



7 ^ 



l - exp(-z) 1 + E rr 



n=l n! 



(A3) 



where the total luminosity of the bulge component is normalized, L B = 7 .22nI eB R^ B = 1, 
and 2 = 7.67(i?/i? e s) 1 / 4 . The exponential law curve-of-growth is: 

Fd(R) = L D [1 - (1 + exp(- J R//i)], (A4) 

where the total luminosity of the disk component is set to the disk-to-bulge ratio, L B = 
2nl h 2 = (D/B), if a two-component model is considered, or normalized, L D = 1, if an 
exponential only model is used (in this case L B — 0). 

Both laws are seeing convolved with a 7 = 1.6 psf, following the technique described by 
Saglia et al. (1993). The Fourier transforms of the 7 psfs are given by: 

pJk)= 2irRJ (kR)pJR)dR = exp [-(&&) 7 ], (A5) 
J 

where J (kR) = J 2?r exp (ikR cos 9) d9 is the zero-order Bessel function. The 7 = 1.6 psf 
reproduces well the stellar profiles measured with the telescopes and setups used in Paper 
III (see Saglia et al. 1993). 

A grid of seeing-convolved models is obtained for 100 values of the V / R eB and F/h ratios, 
ranging from 0.01 to 1 with linear increment of 0.01. Here T is the FWHM of the seeing 
profile. For each of these values, the seeing-convolved luminosity profiles I B (R/T,T/R eB ) 
and I%(R/r,F/h) and curves of growth F^(R/T,T/R eB ) and F£(R/T, V/h) for both the 
bulge and the disk component are tabulated for < R/Y < 50 on a logarithmic radial 
grid (plus R — 0) with dlnR/T = 0.230258 and 31 points. A cubic spline interpolation on 
\n(R/T) and a linear interpolation on T/R eB or T/h are used to determine the profile at a 
given radial distance R and with given values for R eB) h and T. When R/Y < 0.05 a log-log 
extrapolation is used. When R/Y > 50 the correction computed for _R/T=50 is applied. If 
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T/R eB < 0.01 (T/h < 0.01) the correction computed for T/R eB = 0.01 (T/h = 0.01) is used. 
When T/R eB > 1 (T/h > 1) the correction computed for T/R eB — 1 (T/h = 1) is applied. 
The resulting numerical errors in the seeing convolved bulge and disk models are negligible 
(« 1%). 

The luminosity profile Jb+d — /b + jo fitted to the data takes into account the effect of 
the finite pixel size of the observed profiles. These are computed as the azimuthally averaged 
flux in the annulus of radius R and of half pixel width. This procedure is reproduced by the 
following equations: 



f B (R, R e B, T, S) — — 



1 



R + S/2 T 



R 



eB 



-F. 



R-S/2 T 



ReB 



f D (R,h,T,S) = jlFE 



R + S/2 T 
f 'h 



-F, 



R-S/2 T 
f 'h 



(A6) 



(A7) 



where A = ir[(R + S/2) 2 — (R — S/2) 2 ] is the area of the annulus and S is the scale or pixel 
size in arcsec. Eqs. [AI] and |A7] are valid for R > S/2. If R < S/2 (i.e., the central value at 
R = 0), then: 



fs(R, ReB, T, S) 

f D (R,h,T,S)-- 





'R+S/2 r 


b B 


r ' ReB \ 



ir(R + S/2f 



R+S/2 r 

r ' h 



ti(R + S/2) 2 ' 



(A8) 
(A9) 



Similar tables of seeing convolved profiles were also constructed for the ^ = 12 
model and for the smoothed i? 1//4 used in Saglia et al. (1993). The luminosity profile of the 
^ = 12 /oo model is more centrally peaked than the _R 1//4 and declines less rapidly than the 
R 1 / 4 law at large radii. The smoothed i? 1//4 model is less centrally concentrated than the 
R 1 / 4 law. Both profiles have been used to test our fitting algorithm (see §|3.5|). 
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